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Do Solar system tests permit higher dimensional general relativity? 
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ABSTRACT 

We perform a survey whether higher dimensional Schwarzschild space-time is compatible 
with some of the solar system phenomena. As a test we examine five well known solar system 
effects, viz., (1) Perihelion shift, (2) Bending of light, (3) Gravitational redshift, (4) Gravita- 
tional time delay and (5) Motion of test particle in the framework of general relativity with 
higher dimensions. It is shown that the results related to all these physical phenomena are 
mostly incompatible with the higher dimensional version of general relativity except that of 
Motion of test particle. We compare all these results with the available data in the literature. 

Key words: gravitation - Solar system: general - celestial mechanics. 



1 INTRODUCTION 

Einstein's general relativity (GR) received first instant success due 
to the observational confirmation of two solar system effects, firstly, 
contribution to perihelion shift of 43 arcsec per century as curvature 
effect of space-time, and secondly, total solar eclipse in the year 
1919 which admits the relativistic value 1.75 arcsec as obtained by 
Einstein which is also due to the effect of curvature. Probably these 
observationa l boldness and the sublime structure of GR inspired 
iBornl dl962h to state that, "The theory appeared to me then, and 
it still does, the greatest feat of human thinking about nature, the 
most amazing combination of philosophical penetration, physical 
intuition, and mathematical skill. It appealed to me like a great 
work of art ..." . 

The above mentioned two triumph of GR is obviously 
based on its usual four-dimensional structure of space-time. This 
prompted people to start thinking of the multidimensional structure 
of GR. However, the extension of GR by the inclusion of dimen- 
sions beyond four were initiated by investigators mainly in connec- 
tion to the studies of early Universe. It is commonly believed that 
the four-dimensional present space-time is the compactified form of 
manifold with higher dimensions (HD). This self-compactification 
of multidimens i ons have been t hought of by several researchers 
dSchwarzH l985: IWeinbergl 1 198fj) in the area of grand unification 
theory as well as in superstring theory. 

In the Kaluza-Klein gravitational theory with higher dimen- 
sions, therefore, it is a common practice to show that extra di- 
mensions are reducible to lower one, specially in four-dimension 
which was associated with some physical processe s. Interestingly , 
mass have been considered as the fifth dimension dWessodfl983l : 
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199dJPonce de Leonl 2003 ) in the case of five-dimensional Kaluza- 
Klein theory. iFukuil ( 119870 suggested that expansion of the Uni- 
verse follows by the percolation of radiati on into 4-dimensiona l 
space-time from the fifth dimensional mass. lPonce de Leonl ( 2003) 
argued that the rest mass of a particle, perceived by an observer 
in four-dimension, varies as a result of the five-dimensional mo- 
tion along the extra direction and in the presence of elctromagnetic 
field is tot ally of elctromagnetic origin which has confirmed by 
IrHvI |2006). On the other hand, it have been show n by many in- 
vestigators dlshiharall 19841 : iGegenberg & Das. 19851) that within the 
Kaluza-Klein inflationary scenario of HD a contraction of the in- 
ternal space causes the inflation of the usual space. There are cases 
in FRW cosmologies where th e extra dimensions contract as a re- 
sult of cosmological evolution dlba'nez & Verdaguej [l986). In the 
solution to the vacuum field equations of GR in 4 + 1 dimensions 
IChodos & Detweilerl fl98ch have shown that it leads to a cosmol- 
ogy which at the present epoch has 3+1 observable dimensions in 
which the Einstein-Maxwell equations are obeyed. 

Under these theoretical background, therefore, now-a-days 
people have started to think of the higher dimensional influence 
on GR, more precisely, whether within the framework of higher 
dimensional GR the same type of solar system tests would yield 
the same results. Actually, it has two-fold intentions: firstly, if 
the results are positive then the higher dimensional version of 
GR will prove itself as an extended viable theory of gravitation, 
and secondly, if negative then there is no need of higher dimen- 
sional GR at al l . Mot ivated by this, therefore, in a recent work 
iLiu & Overduinl (2000) argued that to test the theory involving the 
motion of test particles in the field of a static spherically-symmetric 
mass like the Sun or the Earth would be most straightforward. 
Kagra manova. Kunz & Lammerzahll d2006h have investigated So- 
lar system effects in Schwarzschild-Sitter space-time and estimated 
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the valu es for the cos mological parameter A. In a similar line of 
thinking llorid d2005d lbl) attempted to investigate secular increase 
of the Astronomical Unit, perihelion precessions and planetary 
motions as tests of the Dvali-Gabadadze-Porrati multidimensional 
braneworld scenario. 

In this co nnection it is to be noted here that in- 
vestigations by | Liu & Overduinl d2000h. along with those of 
llim. Overduin & Wessonl d 19951) and iKalligas. Wesson & E veritt 
dl995h . are limited to five-dimensional soliton-like space-time only. 
Therefore, our present attempt is to study more general cases un- 
der a spherically symmetric Schwarzschild-like space-time with TV 
number of dimensions where N = D + 2 such that D ^ 2. In this 
context we discuss the following five cases involved in the solar 
system experiments to examine the viability of GR with HD, viz., 
(1) Perihelion shift (2) Bending of light (3) Gravitational Red-shift 
(4) Gravitational time delay and (5) Motion of test particle. Our 
present studies show that most of these solar system phenomena 
do not allow dimensions beyond 4 indicating a gross failure of GR 
with higher dimensional framework. 



2 MATHEMATICAL FORMULATION 

Let us consider a spherically symmetric metric which represents 
a generalized Schwarz schild space-time with higher dimensions 
dMavers & Perrvll 19861) 



dr 

da = ~f(r)dt + j^-y + r d£l D , 



(1) 



where r is a radial coordinate and / is a function of r only. The line 
element dflo 2 on the unit D-sphere is given by 

dQ,D 2 = d6\ 2 + sin 2 6\d92 2 + sir? 8\sir? 8id8- 2 + ■■• 



+ Y\ sin 2 8 n d8 D 2 (2) 



with Q D = 2[7r (D+1)/2 ]/[r(£> + l)/2]. Also, according to Ein- 
stein equations we can write f(r) — 1 — /i/r ^ 1 with the constant 
of integration fi = 16itGM/Dc 2 Qd- 

Now, in principle, in Lagrangian mechanics the trajectory of 
an object is derived by finding the path which minimizes the action, 
a quantity which is the integral of the Lagrangian over time. So, in 
connection to the solar system problem we would like to adopt the 
higher dimensional Lagrangian which can be written as 

. 2 

L = T - V = -ft 2 + j+ r 2 6i 2 + r 2 sm 2 0i0 2 2 + ... 



+r 2 Y[ sin 2 8„6 D 2 ■ 



(3) 



Here dot over any parameter implies differentiation with respect to 
the affine parameter V. 

Now, if we take a cross-section by keeping fixed 9i = 82 = 
... = Od-i = f , so that 9i = 0,i = 1,2,3, ...,D - 1 then the 
Lagrangian takes the form 



d_ ( dL 

ds \dqi 



dL_ 

dqi 



0. 



(5) 



By assuming fi — E and t 2 8d = p, where E and p are the energy 
and momentum of the particle respectively, such that i — Ej f and 
8b = p/r' 2 and hence with these notations equation (4) becomes 



c2 -2 2 
E r p 



f f r 

which, after simplification, can be written in the following forms 



(6) 



r* = Lf + E 2 -^ 



and 



1 f dr \' 
r* \d8D~) 



Lf 



f 



(7) 



(8) 



id / P P 1" 

Again, by substituting 8d — <f> and r = 1/U in equation (8), one 
can write 



dU 



Lf , E 2 



fU 1 



(9) 



Now, if we write equation (7) in the form 

\dt) E 2 J E 2 r 2 y ' 

then one can easily observe that ^ vanishes at r = ro of the closest 
approach to the sun. This at once yields the relationship between 
momentum and energy of the particle as follows: p 2 /E 2 = r 2 /fo, 
where fo = f(r — ro). Hence, the equation of photon becomes 



(dr\ 2 2 fr 2 
\dt) 1 for 2 ' 



(11) 



Thus, the time required for light to travel from ro to r can be 
expressed as 



t(r, r ) 



dr 



1/2' 



(12) 



3 SOLAR SYSTEM TESTS FOR HIGHER 
DIMENSIONAL GR 

3.1 Perihelion shift 

Following the equation (8), motion of planet in the sun's gravita- 
tional field can be written as 

4^V = 4 + ?-4 as) 



For r = 1 /U, we have 

^E + U = ^D + 1)U D + f-(D-l)U D - 2 . (14) 
dcp p 

The solution to this equation (14) is then given by the following 
cases: 



L = -fi 2 + j+ r 2 8 D 2 



(4) 



with light-like particle photon, L = and for any time-like particle, 
L = 1. 

Therefore, in terms of the generalized coordinates qi and gen- 
eralized velocities qi, the standard Euler-Lagrange equations are 



(i): D = 2 

By the use of successive approximation (taking fi = as zeroth 
approximation) we get the solution to the above equation (14) in 
the form 



U = ~(1 + ecosq 



(15) 
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Figure 1. The trajectory of equation (15) by choosing suitably the parameters. 



where I = p /^-- Obviously, the trajectory of test particle, i.e., 
planet is elliptical (see Fig. 1). 

Substituting this on the right hand side for U, we get 



U = ~[1 + ecos( 



(16) 



with lo = (3GM/c 2 l)(f>. Therefore, time period for the planet 
is T = 2tt(3GM/c 2 1) and the average precession can be ob- 
tained as n = 6irGM/clT — 43.03 arcsec per century (where 
I = 5.53 x 10 12 cm, GM/c 2 = 1.475 x 10 5 cm and one 
century= 415T). This value of Mercury's perihelion preces- 
sion rate is very close to some of the available data which are 
43. 11 ±0.21 and 42.98 arcsec per c entury respect ively as obtained 
by Shapiro. Cou nselman & Kind dl976t) and lliu & Overduinl 

feooch . 

(ii): D = 3 

The equation of motion in this case can be written as 



d 2 U 



2/j, 



w +u{1 -f ] 



The solution is given by 

U = U cos{j3i>) + Uicos(3l3ip) 



(17) 



(18) 



where ip — acf>, U\ = — (2/i/16a)(7o << Uo and 



3 2 



1 — (3fj,/a)U with a = 1 — 2/i/p . Hence, the path 



is no longer elliptical (see Fig. 2). 



(iii): D = 4 

The solution of the equation (14) related to motion of planet in this 
case is 

"3 1 



U = 



+ 



5/j 



ro r [ 
1 
120 
3/j 



;(2cOS 2 (j> ~ 1) 



— — (8cos 4 (f) — 8cos 2 (f> ~ 1) 



+ 



P 2 r 2 



1 1,„ a 
(2cos 

2 3 ( 



(19) 



Again, one can observe that the path is no longer elliptical (see Fig. 
3). 

lEinsteinl dl915l) explained the perihelion motion of mercury 
from the general theory of relativity by accounting the unsolved 
amount of ~ 43 arcsec as due to gravitati on being mediated by the 
curvature of spacetime fNordtvedt 2001). Besides this relativistic 
effect other effects due to classical reasons are shown in the Table 
1. Therefore, from the present investigation it is revealed that the 
HD model of general relativity only admit four-dimensional case 
with a precession 43.03 arcsec per century. 



3.2 Bending of Light 

Now we would like to observe how higher dimensional version of 
general relativity do respond on the effect of light bending. Let us, 
therefore, start with the equation (9) which now reads 



dU 



E 2 



U 2 {1 



(20) 
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Figure 2. The plot U vs. (j> for D = 3. 



Table 1. Sources of the perihelion shift 



Amount 


Cause 




(arcsec per century) 


5025.6 


Precession of equinoxes 


531.4 


Gravitational tugs of the other planets 


0.0254 


Oblateness of the Sun 


42.98 ± 0.04 


Relativistic curvature of spacetime 


5600.0 


Total 


5599.7 


Observed 



The above equation can be written in the suitable form as 



dd> 2 2 



(21) 



Now, we solve the equation by successive approximation, starting 
with the straight line (path without gravitating body) as zeroth ap- 
proximation such that U = cos(f>/Ro where <j> — is the point P 
of nearest approach to the Sun's surface. Ideally, Rq would be the 
solar radius. 

Substituting this on the right hand side of equation (21) for U, 
we get 



The solution of the above equation (22) is then given by for the 
following cases: 

Case I: D = even = 2n 
Let us consider the case when D is even and takes the value 2n. 
For this particular situation the solution to the equation (22) can be 
given as 



U 



cos(j> fi(2n 
Ro 



2 2n R 2n 



+ ... + 



cos2n<f> 
-4n 2 + 1 4 

2n C n ~icos24> 
-2 2 + 1 



I Cicos(2n - 2)0 
-(2n -2) 2 + l 

M (2n + l) 2n C n 
2 2n + 1 R 2n ■ 



(23) 



(i): For n = 1 
In this subcase 

U = 



Ro 



GM 
R 2 e 2 



(2 — cos <j>). 



(24) 



For U = 0, we get cos<f) = -0.4244302380 x 10~ J where the 
values for the constants are taken as follows: c = 2.997925 x 10 s 
m/sec, G = 6.67323 x 10~ n SI Unit, M = 1.98 92 x 10 30 Kg 
and 7?o = 6.95987 x 10 8 m. We plot U vs. <j> (see Fig. 4). 
Here the net deflection of the ray is given by 



A0 = 1.741300716 arcsec. 



(25) 



This result is in agreement with the experimental result of observed 
deflection of light by the Sun (see Table 2). 



(22) 



(ii): n = 2 

Here 
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Figure 3. The plot U vs. <j> for D = 4. 
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Figure 4. The plot U vs. <j> for D = 2. 
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Table 2. Deflection of Starlight During Eclipses 



Table 3. Deflection of Starlight for HD-Models 



Date 



Location 



Deflection (A</>) 
(arcsec) 



29 May 1919 


Sobral 


1.98 ±0.16 




Principe 


1.61 ±0.40 


21 Sep 1922 


Australia 


1.77 ±0.40 






1.42 to 2.16 






1.72 ±0.15 






1.82 ±0.20 


9 May 1929 


Sumatra 


2.24 ±0.10 


19 June 1936 


USSR 


2.73 ±0.31 




Japan 


1.28 to 2.13 


20 May 1947 


Brazil 


2.01 ±0.27 


25 Feb 1952 


Sudan 


1.70 ±0.10 


30 June 1973 


Mauritania 


1.66 ±0.19 



u 



COS(f> (1 4^ 2(1 

Ro 6-/?q ^-Rq 



4m 



(26) 



For U = 0, we get cos(j> ^imaginary. We plot U vs. 4> (see Fig. 5). 

(iii): n — 3 
For the value n = 3, we get 

ostf> 7 a 



U 



R 



i) 



-L(32co/ 



4:8cos 4 (j} + 8cos 2 c 



+ ^(8cos 4 <j 



8cos 2 cj> + 1) + 5{2cos 2 </> - I) 



+ 



-1) 

35(i 
32Kj 



(27) 



Here (7 = yields cos<?(> = -0.2255226204 x 10" 41 so that 
A(f> = -0.948825313 x 10" 2 arcsec. This value of A0 expresses 
the angle of surplus rather t han angle of deficit jDver & Marleaul 
ll995l : lRahaman et al.ll2005ah . We plot U vs. <f> (see Fig. 6). 



Case II: D 
U-- 



odd = 2n ■ 



+ 



1 

(in 



Ro 



22n . 2R 2n-l 



2n - 1 C 1 cos(2n~3)(t) 
+ -(2n-3) 2 ±l + ' 



+ (in 



cos(2n — l)<f> 
-(2ra - l) 2 + 1 

2n_1 C n -i </>sm</> 

2 2„-2_ R 2n-l 2 



(i): n = 2 
In this subcase 



U = 



cos 



2 jit 

Ro + M 



sirup - 7^cos3(j> 



(28) 



(29) 



For (7 = 0, we get cos^) = -0.2255226204 x 10" 41 . We plot (7 
vs. </> (see Fig. 7). Here the net deflection of the ray is given by 

-2 



A<j) = -0.948825313 x 10" " arcsec. 
which is nothing but angle of surplus. 



(30) 



(ii): n — 3 
In this subcase 

cosd> 3u 

u — — — ^ — — 



^ ■ A 5 ^ 1 KA 



(3D 



For U = 0, we get cos</j = -0.2255226204 x 10~ 41 . We plot U 
vs. (see Fig. 8). Here the net deflection of the ray is given by 

-2 



Aij>= -0.948825313 x 10"" „,-<■.■«;; 
which is again angle of surplus. 



(32) 



Dimensions (TV) 


Deflection (A</>) 
(arcsec) 


4 


1.741300716 


5 


-0.948825313 x 10" 2 


6 




7 


-0.948825313 x 10~ 2 


8 


-0.948825313 x 10~ 2 



Thus, we observe that when D = 2, viz., the total dimen- 
sions, N = D + 2 — 4 then only the result does agree with the 
observational data and D > 2, i.e., N > 4 is not compatible with 
solar system (see Table 3). It can be noted, from the Figs. 4-8, that 
all the trajectories of the light rays almost same for all D due to 
the factor cos<j>/ Ro which is the dominating one. The trajectory of 
light will show different graph for large D. 

Historically, it is important to note that on the basis of his 
'corpuscular' t heory including laws of mechanics and gravitation, 
iNewtonl ( fl704h raised the pertinent issue that "Do not Bodies act 
upon Light at a distance, and by their action bend its Rays, and is 
not this action strongest at the least distance?" He calculated the 
amount of bending of light rays for Sun as 2m/ro. For m = 1475 
meters, in the gravitational units, and ro = 6.95 x 10 8 meters 
thi s equals 0.875 arcsec. However, though prediction of bending 
bv lEuisteinl dl91ll) was at first identical to that of Newton but later 
on he (1915) got the angular deflection of light as twice the size 
he predicted earlier which caused due to the general relativistic ef- 
fect of the curved space-time. In 1919 scientific expeditions per- 
formed at Sobral in South America and Principe in West Africa 
by the leadership of Eddington. The reported observational results 
of angular deflections due to the solar eclipse were 1.98 ± 0.16 
and 1.61 ± 0.40 arcsec, respectively. The mean of these two data 
was taken as confirmation of Einstein's prediction of 1.75 arcsec 
(see Table 2). However, the experiments of Eddington and his co- 
workers had only 30 percent accuracy where the resul ts were scat - 
tered between one half and twice the Einstein value Jwiifeooil) . 
An analysis of large amount of Very Long Baseline Interferome- 
try (VLBI) observations has shown that the ratio of the actual ob- 
served deflections to the deflections predicted by general relativit y 
is very close to unity (e.g., 0.9996 ± 0.001 7 jLebach et all] 19951) . 
0.99994 ± 0.00031 jEubanks et al.lll999h . 0.99992 ± 0.00023 
( Shapi ro. Davis, Lebach & Gregory||2004l) ). 



3.3 Gravitational Redshift 

GR predicted that the frequency of the light would be affected 
due to gravitational field and is observable as a shift of spectral 
lines towards the red end of t he spectrum. Pound-R ebka-Snider 
dPound & Rebkal I1959L Il960l : iPound & Sniderl Il964h confirmed 
this effect through their precision test, sometimes known as Har- 
vard Tower Experiment. In their first test they measured the red- 
shift experienced by a 14.4 Kev 7-rays from the decay of Fe 57 for 
a height of 22.5 meter tower and found z = 2.57 ± 0.26 x 10" 15 . 

Now, as usual, gravitational redshift for the solar system can 
be defined as 



A7 

7 



gu(R) 



1/2 



(33) 
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Figure 5. The plot U vs. <j> for D = 4. 
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Figure 6. The plot U vs. $ for D = 6. 
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Figure 7. The plot U vs. </> for D = 3. 



where R is the radius of the sun and R* is the radius of the earth's 
orbit around sun. 

Therefore, in view of the given HD line element (1), the metric 
tensors involved in the above equation (33) reduce to 

1 M "|l/2 



~9u(R*)~ 


1/2 


_ gu(R) 





1 



1 + ■ 



/' 



jD-1 



(34) 

By substituting the expressions of equation (34) in equation (33) 
for the assumption R* >> fi, we get 

A7 



7 



V_ 

2R D ~ 



(35) 



Thus, in the Sun-Earth system we observe that for the usual 4- 
dimensional case (D = 2), gravitational redshift becomes z ~ 
2.12 x 1CT 6 = z 2 (say). Therefore, for D > 2, z < z 2 which in- 
dicates that as dimension increases the redshift gradually decreases 
(see Fig. 9). It can be also observed that redshift gradually increases 
with mass (since z oc /i = IQtxGM / Dc 2 Q,d, when radial distance 
and dimensions remain fixed in equation (35)). Thus, it seems that 
dimension acts as inversely proportional to mass of the gravitating 
body. 

3.4 Gravitational Time Delay 

Gravitational ti me delay, also k nown as Shapiro time delay which 
was reported bv lShapirq j 19641) is basically the effect of radar sig- 
nals passing near a massive object take slightly longer time for a 



round trip as measured by the observer than it would be in the ab- 
sence of the object there. To proceed on towards the 'Fourth Test of 
General Relativity' let us consider the equation (12) in the form 



t(r,r ) = 

(1-P^r) 
which, after simplification, yields 



dr 



D-l 

r o 



1/2 



(36) 



t(r,r )= / (l-^i)- 1/2 x 



^ ' n 1 1 



( r 2 _ r 2y D -l r D-3 



dr. 



(37) 



Hence, transit time of the light ray from Mercury to Earth can be 
given by 



t = 
+ 

act 

to = 



fri 


r 0\-l/2 
r 2 ' 


' 1-0 


H)\-l/2 
r 2 I 



n 1 I 



( r 2 _ r 2) r D-l r D-3 



1 + 



dr 



+ 



In the absence of the gravitational field (viz., /i = 0) one can 



dr. (38) 



(i-^r 1/2 *+/ r2 (i-Jr 1/2 

r J r 



dr. 



(39) 
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Figure 8. The plot U vs. </> for D = 5. 



Hence, time delay for a round trip is 



(1 _ r °r 1/2 f M + 2V » : —)dr 

( r 2> L r D - 1 (rg-r2)r°-ir °- 31 



At = 2(t - t ) 



- (1 _I| } -v 2 

Let us consider 



( r 2 _ r 2) r D-l r O-3 



I = 



""0 



M 2 



( r 2 _ r 2) r D-l r fl-3 



dr. (40) 



dr. (41) 



The solution to this equation (41) is then given by for the following 

cases: 



At 



4GM 



1+ln 



4rir2 



(44) 



The above expression for radar echo de lay is in accordance with 
the standard literature jWeinberdl2004T) when the Schwarzschild 
space-time is of usual four-dimensional entity and provides an 
amount 240 usee as the maximum excess time delay for the 
Earth-Mercury system. 



(ii): D = 

Here 

3/Lt 

I = sec 

2r \ r , 

Hence, the time delay in this case becomes 



(45) 



(i): D = 2 

In this case the integral in the equation (41) becomes 



I — filn 



,2 _ r 2 





jr — ro 


J r 4 


V r + r 



so that the time delay for a round trip can be given as 



At = — —In 



(n + yVf -rg)(r 2 + yVf - rg) 



+ 



2GM 



ri - r / /' j - 
ri + ro 



r2 + ro 



If, however, ro << ri and ro << ri, then 



(42) 



(43) 



At=^ 
ro 



sec I — I + sec — 

ro / \r 



(46) 



Let us consider that either x = (ri/ro) >> 1 or a; = (ri/ra) » 
1 so that, after neglecting the higher order terms like l/x 3 , l/x 5 , 
... etc. we get 



-l ~ w 1 

sec a; = 

2 x 



so that 



At £S 



4GM 



ro 7r V ri r 2 / 



(47) 



(48) 



(iii): D = 4 
For this case we have 
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I = 



3GM 
2tvc 2 



rro 



- £ + 4 r-ro 
rf ) Vr + r 



(49) 



Therefore, the expression for time delay becomes 

3GM 



At 



c 2 r\ 



n 



T'Z 



r\ - rp 
ri + ro 



+ 



r2 - r 



ri + ro 



(50) 



Thus, from the above case studies one can observe that the max- 
imum time delay will occur when D = 2, i.e., for the usual 4- 
dimensional Schwarzschild space-time. Time delay decreases due 
to increase dimensions. This can be shown easily by assuming 

ro << r\ and ro << r- 2 . 



3.5 Motion of Test Particle 

Let us consider a test particle having mass m which is moving 
in the gravitational field of a D + 2-dimensional spacetime de- 
scribed by the metric ( 1). So, the Hamilton-Jacobi (HJ) equation 
for the test particle is (Chakrabortv 1996; Chakrabortv & Biswas! 
ll99dlRahaman et alj|2005bl) 

ik dS OS 2 ^ 
9 -q— £ + m — (51) 

ox 1 ox K 

where are the classical background field and S is the Hamilton's 
characteristic function. For the metric (1) the explicit form of HJ 



equation (51) now takes the form as 



/ V dt ) J \ dr 



dxi . 



+ ... + 



f dS 



ydxD-i 



\dx 2 ) 



+ m = (52) 



where xi,x 2 , , xb-2 are the independent coordinates on the 

surface of the unit (D — 2) sphere such that 



dQ%,_2 = dx\ + dx\ + + dx 2 D _2 

= del + sin 2 0i_d0l + + 



sin 2 9i sin 2 62 ■ 



. sin 2 6 D - 3 d6 2 D _2 



(53) 



and / is given by f(r) = 1 — /i/r ^ 1 as introduced earlier. 

In order to solve the above partial differential equation (52), 
let us choose the HJ function S as 

S = -Et + Sl(r) +P1.X1 +p 2 .X 2 + +PD-2-X D -2(54) 

where E is identified as the energy of the particle and 

pi,P2, ,Pd-2 are the momenta of the particle along different 

axes on the (D — 2) sphere with the resulting momentum of the 

particle, p = y/p 2 +p 2 2 + + p 2 D _ 2 - 

Now, substitution of the ansatz (54) in equation (52) provides 
the following expression for the unknown function Si which is 



E 2 
f 2 



(55) 
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with e = ±1, where the sign changes whenever r passes through a 
zero of the integral (55). 

To determine the trajectory of the particle following HJ 
method, let us consider that ■§§ = constant and 4^- = constant 

' oE dpi 

[i = 1, 2, , (D — 2)]. Here we have chosen the constants to 

be zero without any loss of generality. 

Therefore, based on the above assumptions one obtains the 
following two integrals 



t 



E 



/ E 2 m 2 p' 2 



' Pi > 



dr, 



i dr. 



(56) 



(57) 



r 2f 



The radial velocity of the particle is then given, from equation (56), 
by 



dr 
~dt 



£2 

7* 



V" 

-prj 



E 



(58) 



The turning points of the trajectory can be characterized by ^- — 
and as a consequence the potential curve becomes 

2 1 1/2 



E 
m 



V7 



V(r) 



(59) 



so that one can write the effective potential, V(r), in the form 

„2 



(60) 

Now, in a stationary system of energy E, the effective potential V 
must have an extremal value. Therefore, the condition to be im- 
posed on the value of r for which energy attains its extremal one 
can be given by 2J- = so that 

u(D - l)mV -p 2 fj,(D+ 1) = 0. (61) 



2p 2 r D 



It has at least one positive root the last term being negative (for 
D > 3). Thus, particles can be trapped by gravitational field 
of higher dimensional Schwarzschild space-time and hence the 
gravitational field is attractive in nature. For D — 2 and D — 3, 
we have some restrictions to get bound orbit as p 2 > 3u 2 m 2 and 
p 2 > um 2 respectively. The plot V 2 vs. r for D = 2 and D — 3 
have been provided in Figs. 10 and 1 1 respectively. 



4 CONCLUSIONS 

Our analytically performed solar system tests for GR with HD can 
be summarized as follows - 

1. Perihelion shift: In 4D our result exactly coincides with that 
of Einstein's predicted value with an elliptical path followed by the 
planet Mercury. As we go increase on dimensions the paths rapidly 
become irregular in shapes and hence HD do not work at all. 

2. Bending light: Here also our theoretical result is in good 
agreement with the experimental result 1.741300716 arcsec which 
become enormously different with an angle of surplus value 
-0.00948825313 arcsec D > 2. 

3. Gravitational redshift: We observe that in the 4-dimensional 
case gravitational redshift becomes z ~ 2.12 x 10 -6 in the Sun- 
Earth system. However, for D > 2 redshift gradually decreases 



with the increase of dimensions such that dimension acts as in- 
versely proportional to mass of the gravitating body. It can also be 
observed that for constant radial distant and dimension the redshift 
gradually increases with the mass of the planets. 

4. Gravitational time delay: It is seen from the present inves- 
tigation that radar echo delay is as usual in the case of 4D and 
decreases with increase of dimensions. 

5. Motion of a test particle: Here the observation is that parti- 
cles can be trapped by gravitational field of higher dimensions and 
hence the gravitational field is attractive in nature (with the restric- 
tions to get bound orbit as p 2 > 3/i 2 m 2 and p 2 > um 2 for D = 2 
and D = 3 respectively). Therefore, this is the only case under our 
study which is fairly compatible with the HD version of GR. 

In a nutshell, our overall obse rvation regarding HP re alm of 
GR is, in general, similar to that of iLiu & Overduir] |2000) which 
is as follows: "... the existence of small but potentially measurable 
departures from the standard AD Einstein predictions". However, 
in some of our HD cases invoke the word 'drastic' in place of 
'small' one! 
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Figure 10. The plot V 2 vs. r for D = 2. 




Figure 11. The plot V 2 vs. r for D = 3. 
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